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Progress in Computing Fuel-Optimal Orbit Transfers 
in Large Numbers of Burns 


C.-H. Chuang and Troy D. Goodson 
Georgia Institute of Technology, Atlanta, Georgia 30332 


Abstract 

This report describes the current state of development of methods for calculating 
optimal orbital transfers with large numbers of bums. Reported on first is the homotopy- 
motivated and so-called Direction Correction method. So far, this method has been 
partially tested with one solver, the final step has yet to be implemented. Second is the 
Patched Transfer method. This method is rooted in some simplifying approximations 
made on the original optimal control problem. The transfer is broken up into single-bum 
segments, each single-bum solved as a predictor step and the whole problem then solved 
with a corrector step. 

I. INTRODUCTION 

Electric propulsion, with its high specific impulse, promises very low fuel 
consumption but it produces less thrust than its counterparts. If one wants to use electric 
propulsion, one needs to be prepared to tolerate the long transfer times that will likely be 
incurred. The greater time spent thrusting must be spent wisely if fuel savings are to be 
realized. Furthermore, the effects of Earth’s oblateness and atmospheric drag become 
more significant on the orbits of long transfer times. 

To spend the thrusting time wisely, form an optimal control problem to maximize the 
mass at the end of the transfer. This, therefore, is the cost function 

J — m(tj) (1) 

subject to the boundary conditions 

v(r(0),v(0),r{r / ) f v(r / )) = 0 (2) 

and the state dynamics 
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r = v 
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v = 
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where ej is the thrust direction, a unit vector, and the thrust magnitude, 7, is limited 
between zero and some maximum value Tmax, fj. is the gravitational constant, g Q is the 
gravitational acceleration at sea-level, and I sp is the specific impulse of the motor. 
Sometimes g(Jsp is referred to as the exit velocity of the motor. 

This results in the well-known bang-bang optimal control problem, discussed in detail 
by Lawden 1 . However, where the boundary conditions are often designed for the 
rendezvous problem, herein the boundary conditions are designed such that the initial and 
final points lie on the desired orbits without specifying the position, or true anomaly, on 
either orbit. 

As found using the Euler-Lagrange necessary conditions, the optimal thrust 
direction for this problem is 



( 6 ) 


where X v is found from the following differential equations 



(7) 

( 8 ) 

(9) 


The optimal thrust magnitude for this problem is a bang-bang solution. Polarity for the 
on-off control is determined by applying the following switching law, Eqn. (10), to the 
switching function, Eqn. (1 1). 


H s > 0 , 7 = 7 ^ 

H s < 0, 7 = 0 
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Solutions of this problem with long transfer times and, therefore, large numbers of 
bums are desired. There are many methods that have been successfully used to compute 
rt-bum transfers, where n is anywhere from 1 to about 6. However, fewer methods 
successfully compute transfers for larger values of n. Methods for the former attempt to 
solve the optimal control problem either directly, indirectly, or as a hybrid of the two. In 
this report, assume that a mostly indirect method, such as BOUNDSCO 2 or MBCM 3 or 
that of Brusch 4 , et. al, or of Redding 5 is being used. 

One idea to obtain interesting solutions is to first compute some n-bum transfer, 
where n may be less than the number of bums desired. Using this as a starting point, 
increase the allowed transfer time and compute the new transfer. It is expected that the 
new transfer is relatively similar to the starting transfer. New transfers are then 
successfully produced this way until the desired transfer is reached. This homotopy 
method seems to work well as long as the number of bums performed in the transfer does 
not need to increase so that optimality is satisfied. For example, in many cases 
BOUNDSCO is unable to find a three bum solution when the two bum solution to an 
almost identical problem is given as the initial guess. Introduced in this report, the 
Direction Correction Method is an attempt to alleviate this difficulty. Its purpose is to 
find an n + 1 bum solution to an orbit transfer problem with allowed transfer time iy+ dis- 
using an n bum solution to the same problem but with allowed transfer time jy . 

Another idea is to patch together a set of n-bum transfers, where n is a small integer, 
perhaps unity, to produce an m-bum transfer, where m is the desired number of bums. 
This method requires that the sequence of transfer orbits be either guessed and iterated 
upon for optimality, or simply prespecified. From the theory of optimal control, this 
patched solution will be a suboptimal solution. This idea will be referred to herein as the 
Patched Transfer method. 

More than likely, once an optimal transfer has been computed, interest will shift to 
developing a guidance law. Possible analytical solutions found from consideration of the 
patched transfer method for the one bum solution of two very close orbits may give a 
simple guidance law. 
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II. direction Correction method progress 


The first idea, referred to herein as the Direction Correction method, is based on 
the common homotopy strategy. The Direction Correction is designed to aid a homotopy 
strategy in calculating successive optimal transfers. In particular, the difficulty arises 
when the desired transfer has one more bum arc than the current computed transfer. 

The method is attractive because it only requires the solution of a relatively small 
set of nonlinear equations. These equations are of the following form 


dD 

dz 


dC 

dz 

i'f) 


m 


<S>(0,t a +dt a )5z(t a +dt a ) = 0 


®(tb + dt b , t f )f (Sz{t a +dt a )) = ~ z(tf)dt f 


dz\ 


<'/) 


(12a) 

(12b) 


For reasons given in a previous report and a paper submitted to the 1994 AIAA 
Guidance, Navigation, and Control conference 6 , both equations are evaluated at time t a + 
dt a . The first equation propagates a guess made for this instant in time back to the initial 
time, using it to check a condition on the boundary conditions at the initial time, denoted 
C(z) where z is the state vector. The second equation is a similar situation, except that it 
is applied to the boundary conditions at the final time, denoted D(z). The function f(z) 
takes into account the fact that the number of bums in the desired transfer, z(r) + 8z(r), is 
one greater than in the computed transfer, z(r). 

The solution information can easily be put into a form useful for a variety of 
numerical methods. For example, the change &(0) can be propagated through the 
transition matrix to calculate the changes at each node point for a multiple point shooting 
method. This method is still under development but shows promise as relatively simple 
way of getting to the n+1 bum solution. 

Using the IMSL routine DNEQBF to solve Eqns. (12a-b), the method has been used 
to predict the correct change, or ‘direction,’ for an example. The algorithm starts with 
information from a given transfer. Then, it iteratively improves upon an initial guess, 
using DNEQBF. The method has produced an approximate solution to Eqns. (12a-b). 
Comparing this solution to the correct answer, errors are only about 4%. The final step is 
to add the solution to the computed transfer and attempt to converge the desired transfer 
with a solver such as BOUNDSCO. Although this has not been implemented yet, success 
is expected. 
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in. patched transfer method progress 


The second method was inspired in part by the work of others. Zondervan, et. al 
made some simple guidance observations 7 , specifically that in some regions the primer 
vector is relatively constant in a velocity-fixed reference frame. This implies that a 
simple control law is available in some cases. Marec presents a solution to the orbit 
correction problem 8 . This motivated a notion that solutions to linearized and/or 
approximated problems were available. In this spirit a solution was obtained for the 
optimal transfer between two close orbits. This solution has been presented in [6]. 

Most interesting about this transfer was the simplicity of the control. Over this 
short transfer between a circular orbit and a close target orbit, the optimal control of the 
thrust angle was almost linear in time. And, in addition, the control direction was almost 
coincident with the velocity direction. 

To review, a modified optimal control problem is considered. The dynamics for this 
problem are again the equations of orbital motion; however, this time the state is defined 
relative to the initial orbit. Assume that the distance from a reference orbit is small 
compared to the radius of the reference orbit and ignore all terms to the order of ( Sr Ip )2. 
This assumption results in the following dynamics: 


8r = 8v 



u(8r^p) 

3 „5 



T 

m = 


(13a) 

(13b) 

(13c) 


Here, 5c and <5v represent displacement from an osculating orbit or the initial orbit, e T is 
the thrust direction, T is the thrust, m is the mass, /i is the gravitational constant, and p 
represents the initial orbit which satisfies identical dynamics but without the thrust term. 
Writing the Hamiltonian for this, the approximated system, gives 


H = X75\ + XJ 


T Cr+3 H(SLtP) 4& 

m p 5 p 3 _ 




o sp 


(14) 


Evaluating the Euler-Lagrange equations results in the following differential equations 
involving the costates: 
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(15a) 

(15b) 
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The control, e T is 



(15c) 


(16) 


and the control T is bang-bang, governed by the switching function, Hj, as 



m 



Sj sp 


H T >0, T - 
H T < 0, 7 = 0 


(17) 


(18) 


Pleasantly, Eqns. (15) happen to be the differential equations for the costates on a 
coast arc coinciding with the initial orbit or the osculating orbit. The coast arc costates 
have been solved by Lawden and other authors 9 * 10 . However, it is Glandorf’s 
formulation, actually based on work by Pines 11 , that is currently being considered. His 
formulation is in the following form: 


MO 

.-MO. 
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(19) 


Considering the form of the state dynamics, their solution can then be written as 


8r(r)- 

8v(0 
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An analytical expression for the integral has been rather elusive. Currently, work is 
focused on approximating the integral. For example, if the magnitude of the Lagrange 
multiplier is approximated as 


fo ) 
-M«.) 


( 21 ) 


where the function g(t) represents a “curve fit” of sorts, then Eq.(20) becomes 



[P(T)] 1 

0 * 
[I 0]P(T) 

ml 
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Now, the integral only has to be evaluated once for each choice of bum, saving a 
considerable about of computation time. Finally, note that bums are not restricted in 
length, using osculating orbits (much as in Encke’s method) the bum lengths are actually 
rather arbitrary. The only consideration for bum length, then, is the error accumulated by 
approximated functions during integration. 

To formulate a method for computing the transfer, the above discussion hints to a 
bum-by-bum approach. Bums would be guessed by a user via a set of transfer orbits and 
bum times. Each bum would then be approached as a single-bum rendezvous problem. 
This produces a sub-optimal transfer and can be thought of as a predictor step. The 
corrector step would then consist of iterations to make it an optimal transfer; either a 
direct optimization of the transfer orbit elements and bum times or an indirect 
optimization by multiple- shooting. 

IV. CONCLUSIONS 


The development of the Direction Correction method is proceeding rather well. 
At the time of this report we are not prepared to say whether the method will be 
successful. The ideas that it is based upon have been validated individually. It has also 
produced a fair approximation to the solution of a known problem. Further testing of the 
method is required in order to determine just how robust it is; but at this point it seems 
pretty clear that method will work. 

The Patched Transfer Method is very promising. GlandorFs formulation for the 
Lagrange multipliers been checked numerically and a suitable approximation for the 
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Lagrange multiplier magnitude is forthcoming. The next steps are to refine the predictor- 
corrector idea, code the method, and test it. 
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